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Abstract 

The wax (glaucousness) on wheat leaves and stems is mainly controlled by two sets of genes: glaucousness loci (Wl and 
W2) and non-glaucousness loci (Iwl and Iw2). The non-glaucousness (lw) loci act as inhibitors of the glaucousness loci (1/1/). 
High-resolution comparative genetic linkage maps of the wax inhibitors Iwl originating from Triticum dicoccoides, and Iw2 
from Aegilops tauschii were developed by comparative genomics analyses of Brachypodium, sorghum and rice genomic 
sequences corresponding to the syntenic regions of the lw loci in wheat. Eleven Iwl and eight Iw2 linked EST markers were 
developed and mapped to linkage maps on the distal regions of chromosomes 2BS and 2DS, respectively. The Iwl locus 
mapped within a 0.96 cM interval flanked by the BE498358 and CA499581 EST markers that are collinear with 122 kb, 
202 kb, and 466 kb genomic regions in the Brachypodium 5S chromosome, the sorghum 6S chromosome and the rice 4S 
chromosome, respectively. The Iw2 locus was located in a 4.1 to 5.4-cM interval in chromosome 2DS that is flanked by the 
CJ886319 and CJ519831 EST markers, and this region is collinear with a 2.3 cM region spanning the Iwl locus on 
chromosome 2BS. Both Iwl and Iw2 co-segregated with the BF474014 and CJ876545 EST markers, indicating they are most 
likely orthologs on 2BS and 2DS. These high-resolution maps can serve as a framework for chromosome landing, physical 
mapping and map-based cloning of the wax inhibitors in wheat. 
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Introduction 

The outermost wax layer protects plants from many types of 
biotic and abiotic stresses, such as drought, phytophagous insects, 
pathogens, solar radiation, and freezing temperatures [1,2]. One 
of the most important roles of the cuticle is to limit transpiration to 
reduce water loss and this provides a key mechanism for plant 
survival in water-limited environments, such as deserts, high 
mountains, saline-alkali lands, and coastal ecosystems [3,4]. 
Worldwide, bread wheat (Triticum aestivum L.) is one of the most 
important food sources for human beings. The wheat leaf, stem 
and, in some cases, spike surfaces are coated with cuticular waxes 
that confer a glaucousness characteristic [5,6]. Physiological 
studies in wheat by Johnson et al. [7] and Richards et al. [8] 
showed that glaucousness reduces transpiration and increased 
water use efficiency. More recently Zhang et al. demonstrated that 
glaucousness reduced cuticle permeability in the terms of non- 
stomatal water loss and chlorophyll efflux [9]. Bread wheat 
cultivars with non-glaucousness traits exhibit significant yield 
increases with reduced solar radiation losses that enable continued 
photosynthesis during the grain filling period [10], and the trait 
may also provide resistance to aphids [11]. 



Glaucousness and non-glaucousness are parallel variations in 
wheat and its relatives. Classical genetic studies have shown that 
both the glaucousness and the non-glaucousness stem and leaf 
phenotypes are controlled by two sets of loci; the wax production 
genes Wl and W2 and the wax inhibitor genes Iwl and Iw2, 
respectively. The Iwl and Iw2 non-glaucousness loci function as 
inhibitors of the Wl and W2 glaucousness loci, and could also 
inhibit other wax production genes in the wax pathway [5,12,13]. 
Genetic analyses have indicated that the Wl wax production gene 
and the Iwl wax inhibition gene are located on chromosome 2BS 
with a genetic distance of 2 cM [12]. However, W2 and Iw2 arc 
separated on chromosome 2DS where the W2 locus is close to the 
centromere [12-15]. Two loci, Iw3 and Ws, were also reported 
conditioning wax on spikes in wheat. Non-glaucousness locus Iw3 
was mapped on chromosome IBS [16] and the Ws gene on the 
short arm of chromosome IAS is responsible for glaucous spikes 
[17]. In addition to these genes, a major QTL (QW.awwSA) that 
accounts for up to 52 percent of the flag leaf glaucousness variation 
has been detected in a doubled-haploid (DH) population [18]. 

Molecular mapping and cloning of genes controlling epicutic- 
ular wax in wheat is of great interests for understanding 
interactions between none-glaucousness genes (lw) and glaucous- 
ness genes (W), as well as their effects on yield, and biotic and 
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abiotic stresses. The Iwl locus originating in wild emmer is closely 
linked to the Xcdo456 RFLP marker at the end of chromosome 
arm 2BS [19]. Liu et al. found that the Iwl locus is 18.77 cM away 
from the powdery mildew resistance gene MIIW170 on chromo- 
some 2BS [20]. Simmonds et al. also reported that the Iwl (Vir) 
gene conditioning a non-glaucousness phenotype maps to chro- 
mosome 2BS [10]. In a tetraploid wheat background, Yoshiya et 
al., have found that Wl is linked to Iwl Dw , but the relationship 
between Iwl and Iwl Du was not confirmed [21], and in an Ae. 
tauschii F 2 segregating population, the non-glaucous locus Iw2 was 
located on chromosome 2DS [22]. In another report, the 
dominant non-glaucous locus Iw3672 (Iw2) derived from a 
synthetic hexaploid wheat also mapped on 2DS by simple 
sequence repeat (SSR) and expressed sequence tag (EST) markers 
[23]. During development of a wheat genetic linkage map with a 
doubled haploid (DH) population derived from the TA4 152-60 
synthetic hexaploid wheat line and the ND495 common wheat 
line, Chu et al., also located a dominant wax inhibitor Iw2 on 
chromosome 2DS [24]. Compared to studies on the Iw non- 
glaucousness loci, little work has been done to map the W 
glaucousness loci in wheat, aside from Wl, which has been 
mapped on chromosome 2BS [21] and the Ws glaucous spike 
allele that is located at the terminus of chromosome IAS [17]. 

The development of a high-resolution genetic linkage map is 
essential for fine mapping and map-based cloning of genes of 
interest. However, it is a tedious undertaking to develop refined 
genetic maps and to clone genes from wheat due to the huge 
genome size (16,000 Mb), polyploidy, high content of repetitive 
DNA and un-availability of a reference genome sequence. 
However, since several grass species including Brachypodium 
distachyon L., rice (Oryza sativa L.), sorghum [Sorghum bicolor L.), 
and maize (geo. mays L.) have been sequenced, and genome 
information from these species provides important resources for 
comparative genomics approaches and for development of high- 
resolution genetic linkage maps of genes of interest in wheat [25]. 
Successful examples of such approaches have been documented 
during cloning of the wheat vernalization gene Vml [26], the 
earliness per se gene Eps-A m l [27], and the durable leaf rust 
resistance gene Lr34 [28], Recently, the shotgun genome 
sequences of hexaploid wheat Chinese Spring [29], Triticum urartu 
[30] and Aegilops tauschii [31] provide more information for marker 
development to the genes interested in wheat. 

In the present paper we have reported (1) identification and 
genetic analysis of the Iwl and Iw2 wax inhibition genes 
originating from wild emmer and synthetic hexaploid wheat and 
(2) development of high-resolution comparative genetic linkage 
maps of Iwl and Iw2 on chromosomes 2BS and 2DS, respectively. 

Materials and Methods 

Plant materials 

Three mapping populations were selected for mapping of the 
wheat wax inhibition genes Iwl and Iw2. To map the Iwl locus, 
WE74, a non-glaucousness common wheat line derived from 
common wheat (glaucousness) and wild emmer (non-glaucousness) 
was used in crosses with Xuezao, a glaucousness common wheat 
line. These crosses produced a 4949 plant F 2 segregating 
population and each F 2 plant was bagged to harvest seeds for F3 
family genotyping. A 120 line DH population developed from a 
hybrid between the non-glaucousness TA4 152-60 synthetic 
hexaploid wheat line and ND495, a glaucousness common wheat 
line, was used to map the Iw2 locus [24]. The newly developed 
International Triticea Mapping Initiative (ITMI) reference map- 
ping population consisting of 1 161 recombinant inbred lines (RIL) 



also was selected for mapping of the Iw2 locus [32]. The 
glaucousness trait was phenotyped on each F 2 plant, F 3 family, 
RILs, and DH lines in field trials with adult plants. Chromosomal 
arm assignment and bin mapping of markers linked to the wax 
inhibition genes Iwl and Iw>2 were carried out with Chinese Spring 
(CS) and homoeologous group 2 nullisomic-tetrasomics [33], 
ditelosomics [34] and deletion lines [35]. 

PCR and product analysis 

Total genomic DNA was isolated from leaves by use of a cetyl 
trimethylammonium bromide (CTAB) protocol [36]. Non-glau- 
cous and glaucous bulks, assembled with equal amounts of DNA 
from 10 homozygous non-glaucous and 10 homozygous glaucous 
F 2 plants, were used for bulked segregant analysis (BSA) [37]. SSR 
and EST markers located on the short arms of homoeologous 
group 2 chromosomes were chosen for polymorphism screening 
[38,39, http://wheat.pw.usda.gov/ cgi-bin/graingenes/browse. 
cgi? class = marker]. Primers for EST markers were designed from 
EST sequences derived from the public NCBI EST database. The 
primer-designed criteria included a Tm of 50-65°C with no 
greater than a 3°C difference between primer pairs. Primer 
sequences and information about the bin-mapped EST markers 
are available at the GrainGenes database (http:/ /wheat.pw.usda. 
gov). Polymorphic markers between the parental wheat lines, as 
well as glaucous and non-glaucous bulks, were genotyped on each 
F 2 plant, RILs, and DH lines. 

Polymerase chain reaction (PCR) was performed in 1 0 (0.1 
reactions containing 10 mM Tris-HCl, pH 8.3, 50 mil KC1, 
1.5 mM MgCl 2 , 0.2 mM dNTPs, 25 ng of each primer, 50 ng of 
genomic DNA, and 0.75 U of Taq DNA polymerase. Amplifica- 
tion of DNA was conducted at 94°C for 5 min, followed by 40 
cycles at 94°C for 45 s, 50-60°C (depending on specific primers) 
for 45 s, and 72°C for 90 s, and reactions were terminated after a 
final extension at 72°C for 10 min. PCR products were mixed 
with 2 ul of loading buffer (98% formamide, 10 mM EDTA, 
0.25% bromophenol blue, and 0.25% xylene cyanol), separated on 
8% non-denaturing polyacrylamide gels (39 acrylamide : 1 
bisacrylamide), and visualized following silver staining. 

Comparative genomics analysis and EST marker 
development 

Polymorphic bin-mapped EST markers flanking the non- 
glaucousness loci Iwl and Iw2 were used in BLAST searches of 
the Brachypodium, sorghum, and rice genome sequences to find 
orthologous genomic regions. Orthologous gene pairs in the 
corresponding genomic regions between the three species were 
located and used to search homologous wheat ESTs (http://blast. 
ncbi.nlm.nih.gov/Blast.cgi) that were used to design PCR primers 
using Primer5.0 (http://www.genome.wi.mit.edu/ftp/pub/ 
software/primerS.O). Polymorphic EST markers between non- 
glaucous and glaucous parental lines, as well as the bulk 
segregants, were used for genotyping three mapping populations 
to construct high-resolution genetic linkage maps. 

Data analysis and genetic linkage map construction 

Chi-squared (% 2 ) tests for goodness-of-fit were performed to 
estimate deviations of observed data from theoretically expected 
segregation ratios. Linkages between molecular markers and the 
wax inhibition loci were analyzed using Mapmaker 3.0 with a 
LOD score threshold of 3.0 [40]. The genetic linkage map was 
drawn with the software Mapdraw V2.1 [41]. 
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Results 

Genetic analyses of wax inhibitors Iwl and Iw2 in 
hexaploid wheat 

At the adult plant stage, leaves and stems of the common wheat 
line, WE74, and the TA4 152-60 and W7984 synthetic hexaploid 
wheat lines were non-glaucous, whereas the common wheat lines, 
Xuezao, ND495, and Opata M85, were glaucous. Fi plants from 
the Xuezao/WE74, ND495/TA4152-60, and W7984/Opata 
M85 crosses were non-glaucous, suggesting that the wax inhibitor 
genes in WE74, TA4 152-60 and W7984 are dominant. The F 2 
population of Xuezao/WE74 segregated as 3730 non-glaucous- 
ness and 1219 glaucousness, which fits a 3:1 ratio (Table 1). The 
F 2:3 progenies segregated as 1269 homozygous non-glaucousness : 
2461 segregating : 1219 homozygous glaucousness, as expected for 
a single gene segregation ratio of 1:2:1 (Table 1). The ND495/ 
TA4 152-60 DH lines segregated identical as reported by Chu 
et al. [22]. The 1161 RILs of the ITMI population segregated as 
549 non-glaucousness and 612 glaucousness to fit the expected 1:1 
ratio (Table 1). These results indicate that non-glaucousness in the 
WE74, TA4152-60, and W7984 wheat lines is controlled by single 
dominant wax inhibitor gene. Since non-glaucousness originates 
from wild emmer in WE74 and synthetic wheat lines in TA4152- 
60 and W7984, the non-glaucousness loci should be designated 
Iwl in WE74 and Iw2 in TA4152-60 and W7984, respectively. 

Identification of SSR markers linked to Iwl 

Because the Iwl and Iw2 non-glaucousness genes are located on 
chromosomes 2BS and 2DS, respectively [10,12,19-24], SSR 
markers assigned to 2BS and 2DS were used preferentially for 
BSA. Linkage of the Xbarc297 and Xwmc25 polymorphic SSR 
markers to Iwl on chromosome 2BS was confirmed after 
genotyping the F 2 segregating populations of the parental lines 
Xuezao and WE74, as well as the non-glaucous and glaucous 
DNA pools of Xuezao/WE74 (Fig. lb). Two SSR markers 
Xgwm614 and Xgwm210, previously linked to Vir or Iwl on 
chromosome 2BS [10], were not polymorphic between Xuezao 
and WE74, or the non-glaucous and glaucous DNA pools, and 
therefore could not be used for Iwl mapping. 

Physical bin map of Iwl and the linked SSR and EST 
markers 

A set of Chinese Spring homoeologous group 2 nullisomic- 
tetrasomics, ditelosomics, and deletion lines [33-35] were 
employed for physical bin mapping of Iwl. The Xwmc25 and 



Xbarc297 SSR markers were located on chromosome 2BS bin 
0.84-1.00 (Fig. la), indicating that the non-glaucousness locus Iwl 
was mapped on the distal part of 2BS. 

The wheat ESTs physically mapping on chromosome 2BS bin 
0.84-1.00 were screened for polymorphisms between Xuezao and 
WE74, as well as the non-glaucous and glaucous DNA pools. One 
EST marker, BE498358, was found to be linked to Iwl (Fig. lb). 

Comparative analysis of the Iwl genomic region 

The BE498358 EST sequence was used as a query to perform a 
Blast search against the Brachypodium, sorghum, and rice genome 
sequences. This search revealed that BradiSgO 1 1 30, Sb06g01410, 
and Os04g0 136700 are orthologs of BE498358 located on the 5S, 
6S, and 4S chromosomes of Brachypodium, sorghum, and rice 
respectively. Putative genes flanking BradiSgO 1 1 30, Sb06g01410, 
and Os04g0 136700 were annotated and compared to identify 
orthologous gene pairs between the three species. The results 
indicate that a 462 kb genomic region in the Brachypodium 
chromosome 5S from Bradi5g01020 to BradiSgO 1430 is syntenic 
to a 3.9 Mb region from Sb06g001110 to Sb06g002790 on 
sorghum 6S and a 5.6 Mb region from Os04g0 118900 to 
Os04g0185100 on rice 4S (Table 2; Fig. 1). Brachypodium genes 
in the syntenic genomic region were then used to find homologous 
wheat ESTs to design primers and for polymorphism screening of 
the parental lines Xuezao and WE74, and the non-glaucous and 
glaucous DNA pools. Ten polymorphic EST markers were 
developed and used to genotype F 2 individuals and to construct 
a high-resolution genetic linkage map of Iwl (Table 3; Fig. lb). All 
1 1 EST markers mapped to the Iwl genomic region, and 8 of 
these could be used simultaneously to find orthologous gene pairs 
in Brachypodium, sorghum, and rice syntenic genomic regions. The 
EST CJ876545 (which is an ortholog of Brachypodium gene 
Bradi5g01160), and ESTs CD927782 and CA695634 were 
excluded. However, the homologous CJ876545 sequences in 
sorghum and rice were not located on chromosomes 6S and 4S, 
but instead are located on chromosomes 10 (Sbl0g005570) and 6 
(Os06g0 182500), respectively. CD927782 was developed from the 
rice gene Os04g0 136900, but the homologous sequence of this 
gene could not be found in the Brachypodium and sorghum 
genomes. CA695634 was located 5.57 cM away from Iwl, and 
thus was not used for further analysis. The Iwl locus co-segregated 
with the BF474014, CJ876545, and CD927782 EST markers and 
is flanked by BE498358 and CA499581 in a 0.96 cM interval in 
the wheat genome that is collinear with a 122 kb genomic region 
with 10 predicted genes that encompasses BradiSgO 1130 to 



Table 1. Genetic analysis of wax inhibitors Iwl and Iw2 in hexaploid wheat. 





Mapping population 


Non-glaucousness* 


Glaucousness* 


Total 


x 2 


X 0.05 


WE74 (Iwl) 


10 










Xuezao 




10 








Xuezao/WE74 


10 










Xuezao/WE74 F 2 


3730 


1219 


4949 


0.36 


3.84 


Xuezao/WE74 F 3 


1269(A)+2461(H) 


1219(B) 


4949 


1.16 


5.99 


W7984 (/w2) 


27 










Opata M85 




20 








W7984/Opata M85 F, 


26 










W7984/Opata M85 RILs 


549 


612 


1161 


3.41 


3.84 



A, H, B represent homozygous non-glaucousness, heterozygous and homozygous glaucousness, respectively. 
doi:1 0.1 371 /journal.pone.0084691 .t001 
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2ISS-4 



Centromere 




Os04gOII8900 



Os04g0 129500 
Os04g0 129600 
Os04g0 129900 



Os04g0131900 

Os04g0 132300 

Os04g0 132400 2 7 
Os04g0 132500 

Os04g0 134800 
Os04g0 134900 



Os04g0 136500 
Os04g0 136700 
Os04g0 136900 
Os04g0137100 
Os04gOI37400 

Os04gOI50300 

Os04g0 154800 



Os04«OI85500 




If 



Figure 1. Comparative high-resolution genetic linkage maps of the wax inhibitors Iw1 and Iw2 in wheat, (la) Physical bin map of Iwh 

Iwl was mapped to distal bin 2BS3-0.84-1 .00. (1 b) Genetic linkage map of Iw1 on wheat chromosome 2BS with genetic distances in cM shown on the 
left, markers shown on the right. (1c) The orthologous genomic region of Iw1 on Brachypodium chromosome 5 with putative genes on the left. (Id) 
The orthologous genomic region of Iw1 on sorghum chromosome 6 with putative genes on the right. (1 e) The orthologous genomic region of Iw1 on 
rice chromosome 4 with putative genes on the right. (If) Genetic linkage map of Iw2 on wheat chromosome 2DS using the ITMI population. Genetic 
distances in cM are shown on the left, and markers are shown on the right. (1g) DH population genetic linkage map of Iw2 on wheat chromosome 
2DS. The markers are shown on the left, and genetic distances in cM are shown on the right. The lines and green solid boxes indicated orthologous 
gene pairs and the blue solid boxes indicated non-orthologous genes between wheat, Brachypodium, rice and sorghum. 
doi:10.1371/journal.pone.0084691.g001 



Bradi5g01220 in Brachypodium, a 202 kb genomic region 
(Sb06g001290 to Sb06g001410) with 13 predicted genes in 
sorghum, and a 466 kb genomic region in rice (Os04g0 132300 
to Os04g0137100) that has 12 predicted genes (Table 2; Fig. 1). 

Comparative genetic mapping of Iwl and Iw2 

EST markers linked to Iwl were used to genotype the 1161 
ITMI RILs and 120 DH lines of ND495/TA4152-60 to develop a 
high-resolution genetic linkage map of Iw2. Eight EST markers 
linked to Iwl were polymorphic between ND495 and TA4 152-60, 
as well as W7984 and Opata M85, and were used to construct a 
linkage map of Iw2 on wheat chromosome 2DS (Fig. 1). The Iw2 
locus co-segregated with EST markers BF474014 and CJ876545, 
and is flanked by the ESTs CJ886319 and CJ5 19831 within a 4.1 
cM sequence in the ITMI RIL population (Fig. If) and a 5.4 cM 
interval in the ND495/TA4152-60 DH population (Fig. lg) on 
wheat chromosome 2DS which is collinear with a 2.3 cM genomic 
region spanning the Iwl locus on 2BS (Fig. 1). 



Discussion 

The aerial surfaces of most plants are coated by epicuticular 
waxes whose chemical and physical properties have important 
roles in interactions between plants and the environment. In wheat 
and its relatives, almost all species have parallel variations of 
glaucousness and non-glaucousness except for Einkorn (A 
genome), which is non-glaucousness [6,12,14]. Genetic and 
cytological studies indicate that glaucousness is mainly controlled 
by two dominant genes, Wl and W2, that are located on the distal 
of 2BS and proximal of 2DS, respectively; and are thought to be 
homologous [12,13]. However, the glaucousness phenotype 
(controlled by Wl and W2) is inhibited by the non-glaucousness 
Iwl and Iw2 loci located on 2BS and 2DS, respectively [12,14,15]. 
These results indicate that the glaucousness locus {W) itself, and 
interactions between the non-glaucousness (Iw) and glaucousness 
(W) loci are responsible for wax phenotypes in different wheat 
tissues. 

Cloning of wheat genes responsible for glaucousness and non- 
glaucousness will provide useful information about molecular 



PLOS ONE | www.plosone.org 



4 



December 2013 | Volume 8 | Issue 12 | e84691 



Comparative Mapping Iwl and Iw2 in Wheat 



Table 2. Colinearity between Brachypodium, sorgh 
IwZ 


um and rice in the syntenic genomic region of wheat wax inhibitors Iwl and 




Wheat EST 


Brachypodium 


Sorghum 


Rice 


Annotation 


CJ5 19831 


Bradi5g01430 


Sb06g001110 


0s04g01 18900 


Hypothetical protein 


CJ9491 74 


Bradi5g01280 


Sb06g001240 


0s04g01 29500 


Sec24-like transport protein 




Bradi5g01270 


Sb06g001250 


0s04g01 29600 


Sec23/Sec24 zinc finger domain containing protein 




Bradi5g01260 






Tyrosine specific protein phosphatase-like 




Bradi5g01250 






Tyrosine specific protein phosphatase-like 




Bradi5g01240 






Hypothetical protein transferase family protein 


CJ777783 


Bradi5g01230 


Sb06g001270 


Os04g0131900 


UDP-glucose:sterol glucosyltransferase 


CA499581 


Bradi5g01220 


Sb06g001290 


0s04g01 32300 


AAR2 family protein 




Bradi5g01210 


Sb06g001310 


Os04g01 32500 


LRR receptor-like serine/threonine-protein 




Bradi5g01200 






N-acylethanolamine amidohydrolase 




Bradi5g01190 






Unknown protein/serine-type peptidase 


BF474014 


Bradi5g01180 


Sb06g001350 


0s04g01 36700 


CBS domain containing protein 


CJ876545 


Bradi5g01160 






UM, zinc-binding; Ubiquitin interacting motif; DA1 -large seed size 




Bradi5g01150 






Plant lipid transfer protein/protease inhibitor/seed storage 


CD927782 






0s04g01 36900 


Conserved hypothetical protein 




Bradi5g01140 


Sb06g001410 


Os04g0137100 


Pectate lyase 15-like 


BE498358 


Bradi5g01130 


Sb06g001410 


Os04g0137100 


Pectate lyase 15-like 




Bradi5g01120 


Sb06g001420 


0s04g01 37400 


Hypothetical protein 




Bradi5g01110 






NB-ARC domain 




Bradi5g01100 






NB-ARC domain 




Bradi5g01090 






NB-ARC domain 




Bradi5g01080 






NB-ARC domain 




Bradi5g01070 






NB-ARC domain 




Bradi5g01060 


Sb06g001460 


0s04g01 47200 


Hypothetical protein 


BQ841470 


Bradi5g01050 


Sb06g001460 


Os04g01 50300 


Hypothetical protein 




Bradi5g01040 






Hypothetical protein 




Bradi5g01030 


Sb06g001490 


Os04g01 54800 


Hypothetical protein 


CJ886319 


Bradi5g01020 


Sb06g002790 


0s04g01 85500 


Zinc finger family protein 


CD882067 


Bradi5g01010 




Os04g0185100 


CDT1-like protein a, chloroplastic-like 


doi:1 0.1 371 /journal.pone.0084691 .t002 



Table 3. EST markers of the wax inhibitors mapped in the Iwl 
and Iw2 genomic regions. 





Wheat 
EST 


Forward sequence (5'-3'| 


Reverse sequence (5'-3'| 


CJ5 19831 


ATACCAAGCCTACTAAGACACTG 


AAGGCATACTCAACAGAAATCA 


CJ9491 74 


TGCTTGGGAATCTGTAATGC 


GCTAACAAATCTGTGGACCTT 


CJ777783 


GCACTCGAAATGACTGGACA 


CACTGCCTTACACTGCAGGA 


CA499581 


GTCACGCTCCTGCTCCTC 


GCACCATCTTGATCCCTCTG 


BF474014 


CCAGTACCTCGAGTCCCAGA 


CGAAGAGGGCGTCGATCT 


CJ876545 


CAAAATGTGATGTCTGCAAGC 


CGATAAGGCCGTTCATATTTGT 


CD927782 TCAGGCAACCAAAACCCTTA 


CCTTTTCTCCAGCTCAATCG 


BE498358 


CAACTACTTCACCCACCACA 


CACTGTGACCCAGGAGCATC 


BQ841470 TTGTTCCGCCTGTATGATGA 


GGAATCCTCATTGGACGAGA 


CJ886319 


GCCATCGGCGTAGTCTTC 


TGGCTTGAAGCAGTGGAAGT 


CD882067 GCGGCAGAAGCTCATATCAT 


TGGGACAAACTCTAGCAGCA 


CA695634 TTAGAAACGACAGTGCAGGG 


GGTGCAAGTACAGAGGAGCC 



doi:1 0.1 371 /journal.pone.0084691 .t003 



interactions between the W and Iw loci, and the mechanisms 
whereby the waxy phenotypes are regulated. Our development of 
a high-resolution genetic linkage map is a first step towards fine 
mapping and map-based cloning of the glaucousness and non- 
glaucousness loci. However, additional refinements to the linkage 
maps are necessary before we can clone the respective genes and 
understand their relationships. 

Comparative genomics analyses have been applied widely to 
develop high-resolution genetic linkage maps of interesting genes 
in wheat [25,42,43] . Macro-colinearity has been observed between 
wheat homoeologous group 2 chromosomes and Brachypodium 
chromosome 5, rice chromosome 4, and sorghum chromosome 6 
[25,44-46]. Several studies have also revealed high levels of micro- 
colinearity in particular genomic regions between wheat, Ae. tauscii, 
Brachypodium, and rice [20,47-51] even through their synteny is 
often interrupted by inversions, deletions, duplications, and 
rearrangements [44,46,47,49]. 

In this study, we have found that a 3.2 cM genomics region 
spanning the Iwl locus in wheat chromosome 2BS was highly 
syntenic to a 462 kb genomic region on Brachypodium chromosome 
5S, a 3.9 Mb region on sorghum 6S, and a 5.6 Mb region on rice 
chromosome 4S (Fig. 1). However, gene duplications, insertions, 
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and deletions were also observed in the syntenic genomic regions 
between wheat, Brachypodium, rice and sorghum. The Iwl co- 
segregating EST marker GJ876545 (orthologous to Bradi5g01 160) 
is not found in the syntenic genomics regions of rice and sorghum, 
indicating that the Brachypodium gene order can serve as a better 
model than rice and sorghum for developing closely linked 
markers in wheat. However, another Iwl co-segregating EST 
marker, CD927782 (orthologous to Os4g0136900), was not 
located in the syntenic genomics region of Brachypodium and 
sorghum, implying that the rice and sorghum genes can provide 
alternative information for marker development of some wheat 
genes. 

The bread wheat genome consists of three subgenomes (A, B, 
and D) that diverged from a common ancestor about 2.5-4.5 
MYA [52,53]. The three subgenomes are still very closely related 
after hundreds of thousands of years of independent evolution and 
genetic linkage maps and comparative analyses over the past 20 
years have revealed substantial conservation of orthologs among 
the A, B, and D subgenomes [38,39]. Conventional genetic 
analyses have also suggested that the Wl and W2 glaucousness loci 
are duplicated genes and that the Iwl and Iw2 non-glaucousness 
loci are also duplicated [12,13]. Additional molecular mapping 
experiments have revealed that both Iwl and Iw2 are located on 
the distal part of chromosomes 2BS and 2DS, suggesting that they 
also may be orthologs. Low polymorphisms are observed on 
chromosome 2DS compared to chromosome 2BS, and of 1 1 EST 
derived markers mapping in the Iwl genomic region, only 7 are 
located on the Iw2 genetic linkage map (Fig. 1). An F 2 mapping 
population containing 4949 plants was used to narrow Iwl to a 
0.96 cM genomic region flanked by the CA499581 and BE498358 
EST markers and this region contained 10 predicted genes in the 
Brachypodium interval from Bradi5g01220 to Bradi5g01 130, the 
same as the result of Adamski et al. [54]. However, by using 1161 



ITMI RILs, Iw2 could only be narrowed to a 4.1 to 5.4-cM 
genomic region corresponding to 26 predicted genes flanked by 
the CJ949174 and CJ886319 EST markers from Bradi5g01280 to 
Bradi5g01020 in Brachypodium. In common wheat, mapping of 
genes on the D chromosomes are often more difficult. The Iwl 
and Iw2 also present such an example. Comparative genetic 
mapping results indicated that Iwl and Iw2 are located in the 
orthologous genomic regions of chromosomes 2BS and 2DS and 
function as orthologs. The mapping results from Iwl will gready 
help identification oiIw2 genes. The EST markers BF474014 and 
CJ876545 are homologous to Bradi5g01180 and Bradi5g01160 
and their co-segregation with Iwl and Iw2 can serve as starting 
points for chromosome landing, physical mapping and map-based 
cloning of the non-glaucousness genes in wheat. 

Acknowledgments 

The authors are grateful to Professor Andrew Jackson, University of 
Galifornian at Berkeley for improving the manuscript. We are also grateful 
to Dr. Steven S. Xu of USDA— ARS, Northern Crop Science Laboratory, 
Fargo, North Dakota, USA, for providing the ND495/TA4152-60 DH 
lines, Prof. Mark Sorrells, Cornell University, USA, for providing the ITMI 
RILs, and Prof. Bikram S. Gill and Mr. W. John Raupp, Wheat Genetics 
Resource Centre, Kansas State University, USA, for providing the Chinese 
Spring aneuploid and deletion stocks. 

All the mapping data, primers are publicly available in the GrainGenes 
database (http://wheat.pw.usda.gov/GG2/ ggdb.shtml). 

Author Contributions 

Conceived and designed the experiments: ZL QS HP. Performed the 
experiments: HWJQJH XZ SO YL DZ ZW. Analyzed the data: HWJQ, 
JH ZWJX. Contributed reagents/materials/analysis tools: HWJQJH XZ 
SO YL DZ ZW QWJX YC. Wrote the paper: HWJH ZL. 



References 

L. Eigenbrode SD, Espclic KE (1995) Effects of plant cpicuticular lipids on insect 
herbivores. Annu Rev Entomol 40:171-194. 

2. Jenks MA, Ashworth EN (1999) Plant cpicuticular waxes: Function, production, 
and genetics. Horde Rev 23:1-68. 

3. Riederer M, Schreiber L (2001) Protecting against water loss: analysis of the 
barrier properties of plant cuticles. J Exp Bot 52:2023-2032. 

4. SchonhcrrJ (1976) Water permeability of isolated euticular membranes: The 
effect of euticular waxes on diffusion of water. Planta 131:159-164. 

5. Jensen NF, Driscoll CJ (1962) Inheritance of the waxless character in wheat. 
Crop Sci 2:504-505. 

6. Bianchi G, Figini ML (1986) Epicuticular waxes of glaucous and nonglaucous 
durum wheat lines. J Agric Food Chem 34:429-433. 

7. Johnson DA, Richards RA, Turner NC (1983) Yield, water relations, gas 
exchange, and surface reflectances of near-isogenie wheat lines differing in 
glaucousness. Crop Science 23:318—325. 

8. Richards R\, Rawson HM, Johnson DA (1986) Glaucousness in wheat: Its 
development and effect on water-use efficiency gas exchange and photosynthetic 
tissue temperatures. Aust J Plant Physiol 13:465-473. 

9. Zhang Z, Wang W, Li W (20 1 3) Genetic interactions underlying the biosynthesis 
and inhibition of a-diketoncs in wheat and their impact on glaucousness and 
cuticle permeability. PLoS ONE 8(1): c.54129. 

10. SimmondsJR, Fish LJ, Leverington-Waite MA, Wang Y, Howell P, ct al. (2008) 
Mapping of a gene (Vir) for a non-glaucous, viridescent phenotype in bread 
wheat derived from Triticum dicoccoides, and its association with yield variation. 
Euphytica 159:333-341. 

1 1. Lowe HJB, Murphy GJP, Parker ML (1985) Nonglaucousness, a probable aphid 
resistance character of wheat. Ann Appl Biol 106:555—560. 

12. Tsunewaki K (1966) Comparative gene analysis of common wheat and its 
ancestral species. II. Waxiness, growth habit and awncdness. Jap J Bot 19:175— 
254. 

13. Tsunewaki K, Ebana K (1999) Production of ncar-isogenic lines of common 
wheat for glaucousness and genetic basis of this trait clarified by their use. Genes 
Genet Syst 74:33-41. 

14. Tsunewaki K (1962) Monosomie analysis of synthesized hexaploid wheat. 
JapJ Genet 37:155-168. 



15. Driscoll CJ, Jensen NF (1964) Chromosomes associated with waxlcssness, 
awncdness and time of maturity of common wheat. Can J Genet Cytol 6:324- 
333. 

16. DubcovskyJ, Echaide M, Giancola S, Rousset M, Luo MC, et al. (1997) Seed- 
storage-protein loci in RFLP maps of diploid, tetraploid, and hexaploid wheat. 
Thcor Appl Genet 95:1 169-1 180. 

17. Gadaleta A, Giancaspro A, Giovc SL, Zachco S, Mangini G, et al. (2009) 
Genetic and physical mapping of new EST-dcrivcd SSRs on the A and B 
genome chromosomes of wheat. Theoretical and Applied Genetics 118: 1015— 
1025. 

18. Bennett D, Izanloo A, Edwards J, Kuchcl H, Chalmers K, ct al. (2012) 
Identification of novel quantitative trait loci for days to ear emergence and flag 
leaf glaucousness in a bread wheat (Triticum ot'stimm L.) population adapted to 
southern Australian conditions. Theor Appl Genet (2012) 124:697-711. 

19. RongJK, Millet E, Manisterski J, Fcldman M (2000) A new powdery mildew 
resistance gene: introgression from wild emmer into common wheat and RFLP- 
based mapping. Euphytica 115:121—126. 

20. Liu ZJ, Zhu J, Cui Y, Liang Y, Wu HB, ct al. (2012) Identification and 
comparative mapping of a powdery mildew resistance gene derived from wild 
emmer (Triticum lurgiduin var. dicoccoidcs) on chromosome 2BS. Theor Appl Genet 
124:1041-1049. 

21. Yoshiya K, Watanabe N, Kuboyama T (201 1) Genetic mapping of the genes for 
non-glaucous phenotypes in tetraploid wheat. Euphytica 177:293-297. 

22. Watanabe N, Takcsada N, Shibata Y, Ban T (2005) Genetic mapping of the 
genes for glaucous leaf and tough raehis in Aegihps iauschii 1 the D-genome 
progenitor of wheat. Euphytica 144:119-123. 

2.3. Liu Q, Ni ZH, Peng HR, Song W, Liu ZY, ct al. (2007) Molecular mapping of a 
dominant non-glaucousness gene from synthetic hexaploid wheat (Triticum 
acstivum L.). Euphytica 155:71-78. 

24. Chu CG, Xu SS, Friescn TL, Faris JD (2008) Whole genome mapping in a 
wheat doubled haploid population using SSRs and FRAPs and the identification 
of QTL for agronomic traits. Mol Breed 22:251-266. 

25. Devos KM (2005) Updating the "crop circle". Curr Opin Plant Biol 8:155-162. 

26. Yan L, Loukoianov A, Tranquilli G, Helguera M, Fahima T, ct al. (2003) 
Positional cloning of the wheat vernalization gene VRN1. Proc Natl Acad Sci 
USA 100:626.3-6268. 



PLOS ONE | www.plosone.org 



6 



December 2013 | Volume 8 | Issue 12 | e84691 



Comparative Mapping Iw1 and Iw2 in Wheat 



27. Faricclli ME, Valarik M, DubcovskyJ (2010) Control of flowering time and spike 
development in eereals: the earlincss per se Eps-1 region in wheat, rice, and 
Br achy podium. Funct Integr Genomics 10:293-306. 

28. Krattingcr SG, Lagudah ES, Spielmeyer W, Singh RP, Hucrta-Espino J, et al. 
(2009) A putative ABC transporter confers durable resistance to multiple fungal 
pathogens in wheat. Science 323:1360-1363. 

29. Brenchley R, Spannagl M, Pfeifer M, Barker GL, D'Amore R, et al. (2012) 
Analysis of the bread wheat genome using whole-genome shotgun sequencing. 
Nature. 491:705-710. 

30. Ling HQ, Zhao S, Liu U, Wang J, Sun H, et al. (2013) Draft genome of the 
wheat A-genome progenitor Triticum urartu. Nature. 496:87—90. 

31. JiaJ, Zhao S, Kong X, Li Y, Zhao G, et al. (2013) Aegilops tauschii draft genome 
sequence reveals a gene repertoire for wheat adaptation. Nature. 496:91-95. 

32. Sorrells ME, Gustafson JP, Somers D, Chao S, Benscher D, ct al. (2011) 
Reconstruction of the Synthetic W7984 x Opata M85 wheat reference 
population. Genome 54:875—882. 

33. Sears ER (1966) Nullisomie-tetrasomie combinations in hexaploid wheat. In: 
Riley R, Lewis KR Chromosome manipulation and plant genetics. Oliver and 
Boyd, Edinburgh, 29-45. 

34. Sears ER, Sears LMS (1978) The telocentric chromosomes of common wheat. 
In: Ramanujam S Proc 5th Int Wheat Genet Symp. Indian Society of Genetics 
and Plant Breeding, New Delhi, 389-407. 

35. Endo TR & Gill BS (1996) The deletion stocks of common wheat. Journal of 
Heredity 87:295-307. 

36. Saghai-Maroof MA, Soliman KM, Jorgensen RA, Allard RW (1984) Ribosomal 
DNA spacer-length polymorphisms in barley: Mendelian inheritance, chromo- 
somal location, and population dynamics. Proe Natl Acad Sci USA 81:8014- 
8018. 

37. Miehelmore RW, Paran I, Kesseli VR (1991) Identification of markers closely 
linked to disease-resistance genes by bulked segrcgant analysis: a rapid method to 
detect markers in specific genomic regions by using segregating populations. 
Proc Natl Acad Sci USA 88:9828-9832. 

38. Roder MS, Korzun V, Wendehake K, Plaschke J, Tixier MH, ct al. (1998) A 
microsatellite map of wheat. Genetics 149:2007-2023. 

39. Somers DJ, Isaac P, Edwards K (2004) A high density microsatellite consensus 
map for bread wheat [Triticum ae.stivum L.). Theor Appl Genet 109:1105 1114. 

40. Lincoln S, Daly M, Lander E (1992) Constructing genetic maps with 
Mapmaker/EXP3.0 Whitehead Institute Tcehn Rep, 3rd edn. Whitehead 
Institute, Cambridge. 

41. Liu RH, Meng JL (2003) MapDraw: a Microsoft Excel macro for drawing 
genetic linkage maps based on given genetic linkage data. Hercditas (Beijing) 
25:317-321. 



42. Guyot R, Yahiaoui N, Fcuillct C, Keller B (2001) In si/ico comparative analysis 
reveals a mosaic conservation of genes within a novel collinear region in wheat 
chromosome IAS and rice chromosome 5S. Funct Integr Genomics 4:47—58. 

43. Schnurbusch T, Collins NC, Eastwood RF, Sutton T, Jefferies SP, et al. (2007) 
Fine mapping and targeted SNP survey using ricc-wheat gene colinearity in the 
region of the Bol boron toxicity tolerance locus of bread wheat. Theor Appl 
Genet 115:451-461. 

44. Sorrells ME, La Rota M, Bermudez-Kandianis CE, Greene RA, Kantety R, et 
al. (2003) Comparative DNA sequence analysis of wheat and rice genomes. 
Genome Res 13:1818-1827. 

45. The International Brachypodium Initiative (2010) Genome sequencing and 
analysis of the model grass Brachypodium distachyon. Nature 463:763-768. 

46. Luo MC, Gu YQ, You FM, Deal KR, Ma Y, et al. (2013) A 4-gigabase physical 
map unlocks the structure and evolution of the complex genome of Aegilops 
tauschii, the wheat D-genome progenitor. Proe Natl Acad Sci USA. 110:7940- 
7945. 

47. Bossolini E, Wicker T, Knobel PA, Keller B (2007) Comparison of orthologous 
loci from small grass genomes Brachypodium and rice: implications for wheat 
genomics and grass genome annotation. Plant J 49:704—717. 

48. Fans JD, Zhang Z, Fellers JP, Gill BS (2008) Micro-eolinearity between rice, 
Brachypodium, and Triticum monococcum at the wheat domestication locus (X Funct 
Integr Genomics 8:149-164. 

49. Lu HJ, Faris JD (2006) Macro- and micro-colinearity between the genomic 
region of wheat chromosome 5B containing the Tsnl gene and the rice genome. 
Funct Integr Genomics 6:90-103. 

50. Griffiths S, Sharp R, Foote TN, Bertin I, Wanous M, et al. (2006) Molecular 
characterization of Phi as a major chromosome pairing locus in polyploid wheat. 
Nature 439:749-752. 

51. Valarik M, Linkiewicz AM, DubcovskyJ (2006) A mierocolinearity study at the 
earliness per se gene Eps-Aml region reveals an ancient duplication that preceded 
the wheat-rice divergence. Theor Appl Genet 112:945-957. 

52. Huang S, Sinkhachornkit A, Su X, Faris J, Gill B, et al. (2002) Genes encoding 
plastid acetyl-CoA carboxylase and 3-phosphoglvccratc kinase of the Triticum/ 
Aegilops complex and the evolutionary history of polyploid wheat. Proc Natl Acad 
Sci USA 99:8133-8138. 

53. Dvorak J. & Akhunov ED. (2005) Tempos of gene locus deletions and 
duplications and their relationship to recombination rate during diploid and 
polyploid evolution in the Aegilops- Triticum alliance. Genetics 171:323-332. 

54. Adamski NM, Bush MS, Simmonds J, Turner AS, Mugford SG, et al. (2013) 
The inhibitor of wax 1 locus (Iwl) prevents formation of f3- and OH-p-dikctoncs 
in wheat cuticular waxes and maps to a sub-cM interval on chromosome arm 
2BS. Plant J. 74(6):989-1002. 



PLOS ONE | www.plosone.org 



7 



December 2013 | Volume 8 | Issue 12 | e84691 



